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Controllable systems relying on quantum behavior to simulate distinctly quantum models so 
far rely on increasingly challenging classical computing to verify their results. We develop a gen¬ 
eral protocol for conhrming that an arbitrary many-body system, such as a quantum simulator, 
can entangle distant objects. The protocol verifies that distant qubits interacting separately with 
the system can become mutually entangled, and therefore serves as a local test that excitations 
of the system can create non-local quantum correlations. We derive an inequality analogous to 
Bell’s inequality [1, 2] which can only be violated through entanglement between distant sites of the 
many-body system. Although our protocol is applicable to general many-body systems, it requires 
finding system-dependent local operations to violate the inequality. A specific example in quantum 
magnetism is presented. 


Quantum simulators can efficiently model quantum 
systems [3, 4]. However, characterizing and validating 
such devices is in general difficult. Indeed, quantum 
state tomography [5] for even eight qubits has required 
weeks of classical computational processing time[6] (in 
addition to exponentially growing measurement require¬ 
ments) . This issue is also present in process tomography, 
the analogue for quantum channels [7, 8]. Notwithstand¬ 
ing the number of measurements growing exponentially 
in system size, for systems with 10+ constituents reli¬ 
able tomography is expected to break down from system¬ 
atic errors in preparation and measurement [9]. Although 
these resource scaling problems are partially alleviated 
with methods based on compressed sensing[10-14], the 
cost still scales at least linearly with the Hilbert space 
dimension, and thus exponentially in the number of con¬ 
stituents. Added difficulties arise when data is limited to 
experimentally accessible local observables [15]. 

Despite the costs of completely characterizing large 
quantum systems, there do exist scalable tests giving in¬ 
complete - though useful - descriptions of system be¬ 
havior. For example, techniques such as randomized 
benchmarking [16-19] and fidelity estimation [20, 21] re¬ 
quire a number of measurements polynomial in system 
size, while still quantifying useful information such as 
error rates or average gate fidelities. In the case of lo¬ 
cally correlated errors, this is sufficient to guarantee the 
operation of error-corrected quantum computer[22-25]. 
Recent linear optics experiments have likewise partially 
verified implementations of boson sampling, a task which 
currently cannot be carried out classically[26]. Although 
it is impossible to efficiently verify this sampling, by us¬ 
ing the results of Ref. [27] the authors of Ref. [28] were 
able to distinguish the sampled distribution from a uni¬ 
form one. The experiment of Ref. [29] produced similar 
results, while also checking whether the photon statistics 
corresponded to indistinguishable particles. 

This paper presents another incomplete test, in the 
context of quantum simulators of many-body systems 


(MBS) [30, 31]. It derives from a constructive procedure 
to entangle two distant ancilla qubits through local inter¬ 
actions with the many-body system. It starts by prepar¬ 
ing the system in an initial, known state, composed of 
many spatially-fixed sites (see Fig. 1). We apply a lo¬ 
cal perturbation to a single site, conditioned on the state 
of an ancilla qubit in superposition. The simulator then 
propagates the system forward in time, correlating the 
ancilla with other sites of the MBS. At one of these dis¬ 
tant sites, we apply a second perturbation controlled by 
a second ancilla qubit. This interaction is chosen to in¬ 
crease the probability amplitude between the current ex¬ 
cited MBS state and its initial state, but only does so if 
both qubits are in the same control state. Qualitatively, 
excitations induced by the first qubit are conditionally 
removed by the second. Although both ancillas are now 
correlated, their correlation with the MBS prevent them 
from being entangled. We therefore measure the MBS 
and post-select for it being in its original state. This both 
disentangles the MBS from the ancillas and increases the 
probability that they are in the same eigenstate. Finally, 
we can directly verify that the qubits are entangled (e.g., 
through state tomography). 

As a toy model, we consider a one dimensional trans¬ 
verse field Ising Hamiltonian, 

= ( 1 ) 

i i 

We write the evolution for a time r as 1^ = ex.p{—irH). 
In the B/J 1 limit, the eigenstates of H are well 
approximated by eigenstates of the transverse field term, 
Hj = —J cF^z^ We therefore assume an initial 

state of the form 

1 ^) = | 0 ) 0 | 0 ) 0 ... 0 | 0 ) , ( 2 ) 

where |0) = |0), which minimizes the energy up 

to corrections of order B • {B/J)^. The lowest lying 
excitations of Hj are described by domain walls. For 
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FIG. 1 : a) and b): Prototype entangling procedure. Ancilla 
qubit A prepared as |+) = (| 0 ) + |l))/\/2 applies controlled 
unitary Gi to site 1 of the initial MBS state \'^)i2e^ producing 
a conditional excited state (|0)^(8)|V^)^2£; + |l)^^i \'^)i2e)/ 
Time evolution Vr propagates the MBS, spreading the exci¬ 
tation across the system. Ancilla qubit B (also in superpo¬ 
sition) then applies controlled unitary U 2 to site 2, removing 
the excitation (conditional on its initial state). Finally, a post¬ 
selection conditioned on projector P is carried out, confirming 
the MBS has returned to its original state. This projects the 
ancillas into a (possibly entangled) correlated state, c) An 
analogous protocol using only one ancilla qubit. After post¬ 
selection, measurements of the single ancilla are sufficient to 
determine whether the previous protocol would have entan¬ 
gled the two ancillas. 

our protocol to produce entanglement, we use excitations 
that are indistinguishable from \g) outside a finite region. 
We therefore consider the evolution of the next lowest 
excitations of Hj^ involving pairs of domain walls of the 
form, 

Kj) = n \9) ’ (3) 

i<k<j 

where 1 < i < j < n. Evolution under H = Hj Hb 
disperses these states across the chain[32]. 

We now consider how to use the MBS to entangle two 
distant qubits. First, an ancilla qubit (labeled A) in 
state 1+) = :^(|0) + |1)) rotates the second} spin of the 

chain through controlled-NOT gate, Ua ,2 = |0)(0U ^ 1 + 

(<2\ 

|l)(l|y^ G) . This produces the excitation |e 2 , 2 ) in the 
chain, conditioned on the state of the ancilla qubit, 

Ua ,2 1+)^ 0 Is() = T (|0)^ 0 \g) + |l)^ 0 |e 2 , 2 )) . (4) 

At this point, the ancilla state is correlated with the 
MBS, though only at spin 2 of the chain. 

Continuing with the entangling procedure, we propa¬ 
gate under W for a time r ^ N/B in order to spread 
the correlation with the ancilla across the spin chain as 


^ Since the first spin only interacts with one neighbor, exciting it 
would create a single domain wall instead of two. 


a superposition of excitations. After applying the \g) 
component of the MBS acquires only an irrelevant global 
phase 6 >, while the excited state |e 2 , 2 ) disperses over the 
subspace spanned by {\eij)}: 

VrUA,2 |+)^®|5) = A |0)^ 0 \g) + |1)^ 0 K • 

( 5 ) 

To understand this process, we consider the projection of 
Hb to the linear span of states an approximation 

valid when BjJ [33]. We note that R b only couples 
between adjacent domain walls, 

(G,iI Hb \ek,i) = 5 (d) 

so the evolution of 1 62 , 2 ) under H is equivalent to that of 
a continuous-time quantum walk in two dimensions [34]. 
The distinct behaviors of \g) and |e 2 , 2 ) under W will allow 
us to distinguish the states of the ancilla through a local 
operation at spin N — 1. 

After evolving for a time r ^ N/B^ the excitation 
1 ^ 2 , 2 ) has a probability ^ 1 /A^^ of being localized in 
state |eAr-i,Ar-i)- We verify this numerically by using 
the quantum walk analogy above. As seen in Fig. 2 a, the 
peak transition probability | (eAr-i,Ar-i| W 1 ^ 2 , 2 ) P scales 
as ^ Key to the success of our protocol, we note 

that the unitary at site — 1 maps the excitation 
\eN-i,N-i) back to the ground state \g). Hence we can 
applying to W |e 2 , 2 ) to give it a non-zero overlap 

with the ground state, 

^ |e2,2) = (ejv-i,iV-i| Vr |e2,2) ^ 0 , ( 7 ) 

where |rp ^ 1/N^. The time r required for the overlap 
|r| to reach its peak scales linearly with N (Fig. 2b), as 
opposed to a time scale ^ N‘^/B observed in diffusive 
propagation [35]. 

We can use the fact that |e 2 , 2 ) can transition to 
|cAr-i,Ar-i) via to entangle ancilla A with a sec¬ 

ond ancilla, B. After time evolution W, we apply a sec¬ 
ond controlled-NOT gate between B and spin N — 

Un-i,b = |0)(0|b O i + |l)(l|s O • (8) 

The resulting state displays correlations between both 
ancillas and the MBS, 

UN-l,BVrUA,2 1 + ^ = 

2 ^ 1^) + \^^)aB ^ 1 ^ 2 , 2 ) + 

|0 1)^5 0 Ifl) + |1 1 )ab ® A^~^^Vr |e2,2)) • 

To motivate the final step, observe that if the ancillas are 
in either state |0 1 )^^ or |1 0 )^^, the MBS is necessarily 
in an excited state. Thus if we measure and post-select 
the MBS to be in the ground state \g), we project the 
ancilla qubits into an entangled superposition of states 
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FIG. 2: Evolution time and ancilla qubit correlation as a 
function of spin chain length N, in the quantum walk ap¬ 
proximation. a.) Inverse plot of the peak values of \r\ = 

I (eiv-i,Ar-i| Vt 1^2,2) I, as a function of chain length N. The 
numerically calculated values (dots) match closely to the lin¬ 
ear fit, with \r\~^ ~ O.OSN + 0 . 55 . b.) A plot of the 
time To taken for the overlap \r\ to reach its peak value. 
The numerically calculated values (dots) match the linear 
ht as tqB ~ 0 . 52 iV + 2.02. c.) A simulation of state oc¬ 
cupations {\eij){eij\) in a quantum walk corresponding to 
N = 32 qubits. The system is initialized in state |e2,2) at 
time T = 0 . The time lapse corresponds to times r = O.Itq 
(top), T = 0.5to (middle), and r = tq, where To corresponds 
to the peak time determined in the previous fit. 

|0 0)y^^ and This post-selection succeeds with 

probability , and produces the entangled state, 

(^1 UN-l,BVrUA,2 \++)aB ^ Id) = 

^==== |0 0)^^ + r |11)^^) . (9) 

^yl + |rP 

To conclude the protocol, a measurement of the ancillas 
would then confirm that the ancillas are entangled. 

This simple approach can be generalized to a proto¬ 
col on generic many-body systems. We consider a prop¬ 
agator for an arbitrary many-body Hamiltonian Vr = 
exp(—i^T) (which may also represent environmental de¬ 
grees of freedom) and an initial prepared ground state 
\i^i2E) (or in the case of a mixed state environment, its 
purification over a larger subspace). The MBS is com¬ 
posed of two local components Si and S2, with E repre¬ 
senting the rest of the system. As in the experimental ex¬ 
ample, we consider two ancilla qubits, interacting locally 
with Si and S2 at different times. In full, the unitary 
evolution is 

Cb{U2) Vr Ca{Ui) 1+)^ |+)b \^12e) , (10) 

where Ca{Ui) = |0){0|a G) 1 + |1)(1U ^ (with Cb{U 2 ) 
defined analogously). Importantly, we assume that the 
local unitaries Ui^U2 individually bring the MBS to an 
excited state, though the combined evolution U2VrUi 
produces a non-zero overlap with the ground state. After 


the second controlled unitary, we make a post-selective 
measurement on the MBS represented by the projection 
operator P. Although this projection can be arbitrary, 
we assume it confirms that the system has returned to 
its ground state. As we argued in the example case, this 
post-selection imparts a known correlation to the ancilla 
qubits. 

The actual test of the MBS derives from verifying that 
the ancilla qubits are entangled. Writing out the qubit 
density matrix in the az basis, we have 

pijrr = I (iulf F iU2y' P {ulY Vr , (11) 

where (O) = (' 0 i 2 £;| O |' 0 i 2 £;) refers to an average over 
the MBS state alone, and p is the probability of making 
the projective measurement P for the state of Equation 
(10). The product (^ 2 )-^ Vr (^ 1 )* (with i^j either 0 or 1) 
represents the unitary evolution of the MBS conditioned 
on the ancillas being in initial state \ij). Hence to gener- 
ically determine whether the qubits are entangled, one 
may carry out full state tomography of the qubits’ den¬ 
sity matrix which can be done through concurrent local 
measurements on the individual qubits. 

The protocol outlined above is the central result of our 
paper. Given the ability to prepare and measure an ini¬ 
tial state, it provides a local test verifying that the MBS 
propagator can generate non-local entanglement. Such 
a result precludes a description in which subsystems are 
locally quantum but all correlations between subsystems 
are essentially classical. In this light, the protocol is akin 
to a Bell’s inequality applied to the system and its dy¬ 
namics as a whole. Like Bell’s inequality, it uses only 
local operations. It also has a ‘loophole’: we require that 
the sites of the MBS are spatially stationary. Otherwise 
a single site could migrate through the MBS and interact 
with both ancilla to entangle them, bypassing the need 
for quantum information to be passed between different 
sites of the MBS. 

Although the procedure we have presented generically 
requires two ancilla qubits to be carried out, under cer¬ 
tain assumptions this requirement may be loosened. In¬ 
deed, although full state tomography on the ancillas is 
required to generically verify entanglement, we note that 
by equation (11) the qubit density matrix is completely 
determined by averages over the MBS alone. This means 
that in certain cases, even in the absence of one of the 
ancilla qubits, it is possible to test whether entanglement 
would have occurred. Such a test derives from the Peres- 
Horodecki criterion[36, 37], which states that the ancillas 
are entangled if and only if the partial transpose matrix 

has a negative eigenvalue. This property is character¬ 
ized by Sylvester’s criterion, which states that a square 
matrix A has no negative eigenvalues if and only if its 
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principal minors are all non-negative[38],^. Hence to con¬ 
firm that the ancilla qubits are entangled, it is sufficient 
to check that a single principal minor of is negative: 
Poi, 01 ^ 10,10 ~ Poi,ioPio,oi < 0. We map this statement to 
an expression on the MBS by using Piji>j> = Pij',i'j 
equation (11): 


v^u^Pu2V^ (ylv^Pvui'^ < (v^Pu2VUi^ 


( 12 ) 

When this inequality holds, the ancilla qubits become en¬ 
tangled under our protocol. It can only be satisfied when 
quantum correlations are propagated between spatially 
distant sites. Importantly, since (O) = {iPi 2 e\ O \i^i 2 E) 
represents averages over only the state \^i 2 e), so it char¬ 
acterizes the many-body system and its evolution alone. 

Although inequality (12) implies the pair of qubits in 
the protocol become entangled, it can actually be mea¬ 
sured using a single ancilla qubit. First, we note that 
the product of terms on the left hand side require no an¬ 
cilla qubits to be measured: for example, the quantity 
{V^U 2 pU 2 V) is simply the probability of measuring the 
MBS in a state corresponding to projector P, after hav¬ 
ing applied the many-body propagator V followed by the 
local unitary U 2 . Contrasting with the left hand terms, 
the right hand side requires an ancilla qubit to measure. 
Preparing the ancilla in state |-h) = ^(|0) + |1)), we 
follow the same unitary protocol as in equation (10), ex¬ 
cept in this case we use the single ancilla as the control 
for both unitaries Ui and P 2 . Written out, this produces 
the state 


\4>) = ^ (|o)^ ® C1^) + |i)^ ® U 2 VU, IV’)) . (13) 


The real and imaginary parts of {V^PU 2 VU 1 ) for the 
many-body state \ 2 p) are then the means of (g) P and 
ay 0 P for the compound state |0), 

(v^PU2VUi'j = {4>\ a^®P\4>)+i {4>\ <7y®P 10) . (14) 

Preparing the state |0) of Equation (13) requires the an¬ 
cilla qubit to interact with both sites of the many-body 
system, but certain cases require only a single site inter¬ 
action to measure {V^PU 2 VU 1 ). This occurs when the 
post-selection projector takes a tensor product form. 


P = P 2 0 Pie 0 Ie' • (15) 


The product PP 2 can also be written in this way, 

PP 2 = (g) Pie (g) . (16) 


^ The principal minors of a square matrix Amn are the determi¬ 
nants det(A(^i’^2, ••• where ••• is the matrix A 

truncated to only rows and columns {si, S2, , s^}. 


where we have written p 2 U 2 in terms of its Hermitian and 
anti-Hermitian parts, and as before we let E' denote the 
(inaccessible) environmental degrees of freedom. Using 
this decomposition, it suffices to prepare the state 

10') = ^ (|0)^ ® C 10) + |l)^ ® VU^ 10)) , (17) 

which requires only a controlled unitary between the an¬ 
cilla qubit and site 1 of the MBS. As before, the right 
hand side of (12) can then be written as a sum of ob¬ 
servables, 

(v^PU2VUi^ = (0'| ®B+-<7y® B_) 10') 

+i (0'| (^ax 0 + cTy 0 5+) 10') (18) 

We note that with identity (11), both of these procedures 
can be generalized to do complete state tomography. 

The ideas we have presented have potential applica¬ 
tions in a variety of existing experimental setups. For 
example, current ion trap experiments have the poten¬ 
tial to simulate spin models displaying long-range prop¬ 
agation of correlations [39-44], making them amenable to 
the single ancilla protocol described above. The models 
studied have ground states that can be both prepared 
and measured through direct fluorescence spectroscopy 
following (if necessary) adiabatic passage. The proto¬ 
col’s ancilla qubit can take the form of either one of the 
ions present in the system or one of the global motional 
modes associated with the ion trap. A key requirement 
in these setups is the ability to individually address single 
ions in the experiment [45-48]. Alternatively, in optical 
lattice many-body simulators [49-53] it is possible to use 
polarization of light as the ancilla qubit. Based on selec¬ 
tion rules arising out of angular momentum conservation, 
the controlled unitary operation of the ancilla would cor¬ 
respond to a polarization-dependent interaction with a 
localized subsystem. Between interactions the light must 
be sent through a delay line (e.g., a Fabry-Perot cavity), 
so that correlations between spatially distant MBS sites 
are given enough time develop[54, 55]. 
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